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Abstract 

The generating function for spanning forests on a lattice is related to the q-state Potts 
model in a certain q —* limit, and extends the analogous notion for spanning trees, or 
dense self-avoiding branched polymers. Recent works have found a combinatorial pertur- 
bative equivalence also with the (quadratic action) 0{n) model in the limit n — > —1, the 
expansion parameter t counting the number of components in the forest. 
We give a random-matrix formulation of this model on the ensemble of degree-fc random 
planar lattices. For fc = 3, a correspondence is found with the Kostov solution of the 
loop-gas problem, which arise as a reformulation of the (logarithmic action) 0{n) model, 
at n = —2. 

Then, we show how to perform an expansion around the t — theory. In the thermody- 
namic limit, at any order in t we have a finite sum of finite-dimensional Cauchy integrals. 
The leading contribution comes from a peculiar class of terms, for which a resummation 
can be performed exactly. 

PACS: 05.50.+q, 75.10.Hk, 02. 10. Ox 
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1 Introduction 

The 0(ri)-invariant cr-model for n = 0, —1, —2 defined on a generic graph G has a very inter- 
esting combinatorial interpretation. 

Already in early 70's, it has been observed that the n-vector model in the limit in which 
n —2 is equivalent to a free fermionic theory [1] . But the quadratic (Gaussian) term is the 
Laplacian on the graph G and the partition function is its determinant, which according to 
KirchhofF matrix-tree theorem, provides, once the zero-mode has been removed, the weight of 
spanning trees of the graph G [2] . 

In 1980 Parisi and Sourlas [3] showed the equivalence of the n-vector model in a limit in 
which n — > 0, which was already known to describe the critical behaviour of polymers [4], with 
a supersymmetric osp(2|2) model in which the loops in Feynman graphs, which vanish in the 
n — > limit, gives zero contribution, because of the cancellation between bosons and fermions 
(independently also McKanc noticed that fermions can be used to cancel the contribution of 
bosonic loops [5]). Also this construction is independent from the choice of the graph G, as 
only exploits the symmetry properties in the target space. The choice for the 0(n)-invariant 



model as a cr-model provides exactly the partition function of self-avoiding walks on the graph 
(see for example [6]). 

More recently, we have shown [7] that the generating function of spanning forests in a graph 
G can be represented as a Grassmann integral of the exponential of a local fermionic weight 
involving a Gaussian term together with a special nearest-neighbour four-fermion interaction 
(see also [8]). Furthermore, the fermionic model possesses a hidden osp(l|2) supersymmetry 
non-linearly realised. In [7] we also discussed briefly how this fermionic model can be mapped, 
at least in perturbation theory, onto an osp(l|2)-invariant tr-model with spins taking values in 
the unit supcrsphere in R-"^!^, or an 0(n)-invariant (j-modcl with spins taking values in the unit 
sphere in (also known as n-vector model), analytically continued to n = — 1. The parameter 
t which appears in the generating function of the forests to count the number of trees in a forest 
is related to the coupling constant in the cr-model (with an important inversion of sign). It is 
remarkable that the same generating function can be obtained by a suitable limit g' ^ of the 
g-state Potts model defined on the same graph G (see for example [9]). 

Very detailed information on these models can be obtained by considering regular graphs, 
and in particular in two dimensions, where methods of Conformal Field Theory and Integrable 
Systems apply. 

The critical exponents of the n- vector model, at least on the range n G [—2,2], can bo 
computed exactly, thanks to a mapping onto the solid-on-solid model [10] of a suitable choice 
of the weights, which, in the high-temperature expansion forbids loop crossing. We shall call 
this variant of the n-vcctor model the Nienhuis m,odel, or Loop-gas model (as it is a 'hard-core 
lattice gas' in which the elementary objects are self-avoiding loops). Also the critical behaviour 
of the Potts model can be analyzed exactly, thanks to the mapping onto an ice-type model, 
which has been constructed both algebrically [11] and combinatorially [12]. This is once more 
mapped to a solid-on-solid model [13]. And the critical limit of the solid-on-solid model is 
recovered by using the Coulomb-gas picture [14]. More precisely, the critical behaviour of this 
n-vector model is given by a conformal field theory (OFT) with central charge 

c{n) = l-—^—— (1.1) 

where the parameter m is related to n by the relation 

n = 2cos— (1.2) 
m 

and to the Coulomb-gas coupling constant by 

90 = 1 + -. (1.3) 

m 

Please, remark that c(— 1) = —3/5 and c(— 2) = —2. This means that the Nienhuis model does 
not describe at n = — 1 the universality class of spaiming forests. Indeed, a direct perturbative 
analysis of the 0{n) cr- model on a square lattice [15] or on a triangular lattice [16] at n = — 1 
shows that the model is asymptotically free for t = 0+. As a result, the ultra-violet fixed point 
is the free theory which describes trees, therefore, aside logarithmic violations, the central 
charge is c = —2. 

Beyond the study of a model on a fixed periodic planar graph, a lot of progress has been 
achieved by considering an ensemble of planar graphs [17,18]. Such a study has both an interest 
per se of combinatorial nature, and a relevance in connection to the original case of regular 
graphs. Indeed, after the work of Knizhnik, Polyakov and Zamolodchikov [19,20], it is now 
understood that, for systems showing conformal invariance at criticality, statistical averages 
in the two ensembles have related critical behaviours, so that informations in one context can 
be inferred from the other, e.g. concerning the critical exponents, the conformal families of 
operators, and their dimensions. [19-23] In many cases, cross-checks of these predictions have 
been performed [23-25]. 
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A deep understanding of KPZ relation is a hard and active field [21,26,27]. However, some 
heuristic reasons can be given at least for the existence of a relation of this kind. In two di- 
mensions, at criticality, scale invariance, combined with (discretized) Euclidean symmetries, is 
promoted to the symmetry described by the full Virasoro Algebra. Analogously, statistical me- 
chanics models on random planar graphs, when reaching simultaneously the largo- volume limit 
and the critical point for the "matter fields" {double scaling limit), show the scale invariance 
pertinent to criticality, combined with the (discretized) invariance under local diffeomorphisms 
(as bare random planar graphs describe a discretization of two-dimensional quantum gravity). 
But a conformal theory in presence of two-dimensional quantum gravity enjoys a symmetry 
corresponding to a SL(2,M) current algebra [19], which is larger than that described by the 
Virasoro algebra. 

This richer structure is in a way at the root of the fact that many results exist for the 
apparently harder counting problem on random planar graphs (which involves a double av- 
erage), and still lacks in the Euclidean case (and indeed provide a hint to critical aspect of 
these quantities, through KPZ). In words more appropriate to the discrete setting, in many 
combinatorial approaches (among which the present paper), the interplay between degrees of 
freedom of the lattice and of matter fields plays a crucial role in simplifying the expressions. 
On the other side, many Euclidean concepts involving a natural notion of distance are harder 
to define cleanly in the random-graph setting (and, if defined through geodesic distance, hard 
to compute). 

On the other side, the generating hmction of statistical configurations over random graphs 
(as well as many other "global" physical observables, such as susceptibilities) can be written as 
the Feynman expansion of a proper action of a zero-dimensional field theory: the replacement 
of real or complex bosonic fields with N x symmetric- or hcrmitcan-matrix fields allows to 
count graphs of genus f) with a weight proportional to N~'^'^, and a large- A limit, achieved via 
steepest descent or continuous approximation of matrix spectra, gives the restriction to planar 
graphs. Such a strategy, started with the seminal works of [17, 18], had a strong development in 
the subsequent thirty years, and now deserves the name of Random Matrix technique (see [28] 
for a recent pedagogical introduction). 

Indeed, after the Ising [29] and Potts [30] models, also the Nienhius model has been 
solved [31,32] on random planar graphs, and deeply studied (see for example [33-35]). In [36], 
as we are interested in the cases n = — 1, — 2 of the Nienhuis model, a combinatorial reformu- 
lation of these problems has been introduced to achieve the random matrix solution with no 
need of an analytical continuation. 

A detailed account of these different research areas when they overlap on geometrical critical 
phenomena is given by Duplantier and Kostov [25]. 

But the model of spanning forests in an ensemble of random graphs escapes the realm of 
exact results. Only the limit of spanning trees, that is f = 0, that we have seen corresponds 
to the Nienhuis model at n — —2, has been studied, and for regular graphs with coordination 
number 3. For spanning forests it is necessary to dispose of informations at finite values of the 
coupling constant, at least in perturbation theory. 

A full discussion on the contents of this paper is postponed to the end of section 2, after 
that some other definitions are introduced. 

2 The model 

Since here on, and according to common use in Random Matrix, we call a "graph" G (or, 
more precisely, a "fatgraph"), what is indeed an orientable 2-dimensional cell complex, i.e. G 
is determined not only by the sets V{G) and E{G), respectively for vertices and edges, but 
also by a consistent choice of the set F{G) for the faces (or, equivalently, for each vertex, the 
outgoing edges have a given cyclic ordering). For this reason, the genus 1^ of G is univocally 
and easily determined, for example via Euler formula {V, E, F and K denote respectively the 
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number of vertices, edges, faces and connected components in the cell complex) 



2l) = 2K + E-V-F. (2.1) 

Consider the ensemble of all connected graphs with vertices of degree k, with a measure 
depending from the genus (and allowing to take a "planar" limit). Most commonly studied 
cases are fc = 3 of fc = 4: we will study the generic case, but with special attention to fc = 3, 
both because it is the ensemble studied in Kostov solution, and therefore this allows for a 
direct comparison of results, and because the generating function Ai(lu) for cubic trees is the 
one with the simplest explicit formula. 

For a connected graph G with V vertices and genus f) we consider the customary (unnor- 
malized) measure for Random Matrix theory 

...(G) = ^^,-A^-", (2.2) 

and the generating function (for connected graphs) is obtained by the Random Matrix technique 
with a single matrix field: 

Zo = E^.^-(G) = l^ln / dMe^K-^'^^+f-^) . (2.3) 

Here the integral is over a set of iV^ real variables, arranged in the Hermitian matrix M, and 
must be intended as the formal Feynman expansion in the parameter g [28]. Indeed, a sketch 
of the technique is the following: the Feynman diagrammatics leads to the generating function 
of all connected graphs, and the traces due to the "matrix of fields" , jointly with the Wick rule 
(MijMgk) — SjiSki, lead to a combinatorics on face-indexing that, via Euler formula, allows to 
count the genus of the graph with the desired factor N~'^^^'^K 

This is a one-matrix theory, i.e. we have only one matrix of fields. We recall here some 
results. First perform the change of variables M — > UAU~^, with U unitary and A diagonal, 
the Jacobian corresponding to the square Vandermonde determinant A^(A) — Yii^j ~ 
Angular degrees of freedom do not appear in the action, and are trivially integrated, similarly 
the ordering of the eigenvalues is irrelevant, and is trivially summed over, and the partition 
function reads 

Zo (X -i^ In / d^A A2(A) exp [a^J] (-^A^^ + f Af) ] . (2.4) 

At this point, many tools allow to extract the asymptotic behaviour in V in the planar limit 
(saddle point, orthogonal polynomials, loop equations. . . ). The result for the leading behaviour 
near to the radius of convergence gc of the series is, for any genus f), 

Zoig,i))-mY^f^yV-'+^+^^'N-'\ j = -l; 7' = ^. (2.5) 

The quantity in (|2.4p is the pure gravity partition function, which should be adopted as the 
appropriate normalization factor when comparing with the analogous expression in the case of 
a matrix theory describing an ensemble of combinatorial structures on the graph (interpreted 
as matter fields coupled to the gravity). In our case, we should deal with Fortuin-Kasteleyn 
random clusters for the analytic continuation in q of the q-colour Potts Model, in a limit g — > 
corresponding to a restriction to forests, with a factor t per tree in the forest. 

So, let T{G) be the set of spanning forests over the graph G. Say for short that F ^ G 
if F £ T{G), and call K{F) the number of its connected components. Given two graphs G, 
H with H C_ G, define G/H the contraction of G by H, i.e. the graph in which vertices 
i and j are identified if an edge {ij) is in E{H) (this is the concept used, for example, in 
deletion/contraction operations for Tutte-Grothendieck invariants). In particular, if G has V 
vertices and E edges, and F ^ G has E' edges, the graph G/F has V — E' vertices and E — E' 
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edges. Furthermore, as any tree is homotopic to a point, G and G/F have the same topology 
(f)(G) = \){G/F)), however they may have dilFerent automorphism groups. 
At fixed graph G, we define the spanning-forest generating function as 

Z(i;G)= ^ Mt,G(i^), (2.6) 

F-<G 

through the (unnormahzed) measure 

„ (p.fK(F) |Aut(G)| 

|Aut(G/F)r (^-^^ 

The symmetry factor, depending on the pair (G, F), has been introduced for later convenience. 
Here we remark that, if we worked with (leg-)labeled graphs, and exponential generating 
function, i.e. 

E (^^^(^)' (2.8) 

G labeled ^ ' 

the rephrasing of the measure above would involve no factors whatsoever. Furthermore, for 
large graphs (which are the dominant class in the thermodynamic limit), these symmetry 
factors are almost surely 1. 

A remarkable exception to this last argument occurs in our exact resummation of section 
|9l where the contributions of certain families of "large-volume" graphs are resummed into the 
evaluation of certain integrals over graphs with size of order 1, for which the inclusion of the 
proper symmetry factor is important. 

We are interested in the macrocanonical average in the random lattice ensemble, mainly in 
the limits N oo {planar limit), and g gdt), the radius of convergence of the resulting series 
(thermodynamic limit). A key fact is that interchanging the two sum operations (over graphs 
G, and over forests F on G) the expression largely simplifies. Indeed, easy manipulations give 

Z{t, g, N)^Y. f'9,N{G)Z{t; G) - ^ ^ M<,,Jv(G)Mt,G(i^) 

G G F^G 

iK[F) 

^^JAut(G/^^)| (2.9) 

^sr- fK{F) \v{F)\ 

F h |Ant(G/f^)| ' 

and in particular, for the planar case TV oo, we have 7V~2i)(G) _^ (5f,(c!) qj and 

planar 

Estimating this function as well as possible, in the neighbourhood of i ~ and g — )■ gdt), that 
is, for graphs in the thermodynamic limit, and near to the critical point of spanning trees, is 
the main goal of this paper. 

In section [3] we give some preliminary results of combinatorial nature. In section |4] we 
describe a standard random-matrix approach to the problem at generic t, g and N. In sections 
[5] and [HI it is shown a correspondence with the 0{n) loop-gas model. 

In sections [THHl the function (|2.10p is shown to admit a diagrammatic perturbative expan- 
sion in the parameter i, i.e. the coefficients Zn{g) of the series 



Z(t,g) = ^rZ„(g) (2.11) 
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can be calculated with a combinatorial technique, each term being a finite sum of finite- 
dimensional integrals. 

The way in which these diagrams are effectively evaluated is described in section [101 while 
sections [H] and [12] describe how to extract the leading behaviour from this sum, in the double 
limit t ^ and g / gc, for t/ \n{gc/g) below a critical threshold. 

While sections [H [5] and [Sj involve concepts of Random Matrix Theory, all other sections 
are purely of combinatorial nature, and use elementary self-contained methods (except for the 
isolated appearence of the non-trivial Levy generalized Central Limit Theorem, in section [Sj 
which furthermore is not mandatory in the logic of the paper, as the results are rederived in 
section [T0|). 



3 Preliminary combinatorial results 

Here we make a short review of results for certain counting problems which will arise in our 
analysis. This section does not make use of any tool from Random Matrix theory, and is fully 
elementary and self-contained. We call H the upper half plane, and D the unit disk. 

Problem 1. Counting the configurations of n non-intersecting arcs in H, with endpoints in 
dm (cfr. fig. [21 top left). 

Such a configuration is called a link pattern (in H). Their number is C„, the n-th Catalan 
number. We adopt the convention Co = 1 for the empty configuration. Call C{q) = ^„ CnO" 
the generating function. A sketch of proof is as follows. For n 7^ 0, one can remove the 
right-most arc. Then, the inner and outer part of the original diagram correspond to smaller 
independent configurations, of sizes n' and n ~ n' ~ 1. This leads to a convolutional relation 
for the coefficients C„, and a polynomial one for the generating function 

C{q) = l + qC{q)\ (3.1) 

the solution matching the regularity condition C_i ~ being 



Cil)-'-^^- (3.2) 



The closed expression for the Catalan numbers is 

ft.^4.r"). («) 

n-\- 1 \ n J 

We say that a tree has degree k if all of its vertices have degree either fc or 1, the latter 
being called leaves (indeed, any tree with at least one edge must have at least two leaves). 
The word cubic is used here as a synonimous of degree 3. We recall that, for a tree to have v^, 
vertices of degree d, the set of u^'s must satisfy the constraint 

Y,{'^-d)vd^2. (3.4) 

d 

Then we will consider the problem 

Problem 2. Counting the configurations of non-intersecting cubic trees, with n vertices of 
degree 3 inW and n -\- 2 leaves in dM (cfr. fig. [3 top right). 

Call Ai^n the number of such trees, and Ai{q) the generating function Ai{q) ~ X]n>o^i^"9"- 
We have Ai,o = 1 for the tree with a single edge and two vertices of degree 1. For n > 0, the 
right-most leaf must be connected to a vertex of degree 3. Then, we can modify the drawing 
in a unique way so that each of the two other terminations can be seen as the right-most leaf 
vertex of the remaining component. Thus we have the formula 

A,{q) = 1 + qA,{q)\ (3.5) 
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Figure 1: On the left, a link pattern on M and the removal procedure which leads 
to the relation HS.lp . On the right, a cubic tree on M, and the removal procedure 
which leads to the relation 13. 5\} . 



i.e. we deal again with Catalan numbers, 

Ai(g)=C(g) = l^^^I^, (3.6) 

and Ai,n = C„. 

We will now consider a generalization to arbitrary degree k: 

Problem 3. Counting the configurations of non-intersecting trees of degree k — h + 2, with n 
vertices of degree fc m H and hn + 2 leaves in dW. 

The reasoning follows as above, with the only difference that now, under removal of the right- 
most leaf and the corresponding internal vertex, we have h + 1 remaining components, thus we 
have 

A^,iq) = l + q{An{q))''^\ (3.7) 
which gives a generalization of Catalan numbers [44] 

Performing a Stirling expansion of (|3.8p for large n, we get 



Conversely, equation (|3.7p for the generating function for general h cannot be written in a 
simple form. A futher exception (besides /i = 1) is the case h ~ 2, for which we have 

, /4a;2\ 3 /I \ , , 

A2 = — sm - arcsmx . (3.10) 

\ 27 / X \3 / 

In Appendix |^ it is shown that the solution to (j3.7p is a generalized hypergeometric function 
h+iFfi, and some properties are studied. Here, we discuss a minimal set of properties which 
are strictly necessary in the forthcoming sections. 
Defining for future convenience the combination 

(3.11) 



7 



we see that the radius of convergence for Ah {q) is 

kl<2NcW. (3.12) 
The series has a finite value also for this value, indeed 

Ah{2'^g,{h))^^, (3.13) 

and, as the series has positive summands, a fortiori \Ah{2'^gcih)e^^) \ < 
Equation (|3.7p has a parametric solution in algebraic form 



gi/'' = z = x{l - x'') 
= (l-x")-i; 



(3.14) 



which is easily checked by direct substitution, and by matching the initial condition A^^ = 1. 
By a simple scaling we have, for z = x{l — gx^), 

Ahigz'^) ^ (l ~ gx'^r' ; z{Ahigz'^) ~ l) ^ gx'^+\ (3.15) 

The problem of counting non-intersecting trees in EI is obviously related to the one of counting 
non-intersecting trees in D (with leaves in dIS>). The solutions essentially do coincide (if one 
point on 9D is marked), except for the fact that, for the case on the disk, it is natural to add 
a symmetry factor for the cyclic permutations of the leaves. Thus we introduce the modified 
quantities 

^''■"=ftfe' AU<?) = ^g"<„. (3.16) 

n>l 

Remark that for A' , differently than for A, wc start summation from n = 1. Again the case 
fc = 3,i.c. /i = l, gives a simple explicit formula 



More generally, the definition p.l6p implies the relation 

z{Ah{gz^)-l) = ^Jz'A'h{gz'^)) , (3.18) 
which, using ^/(z(a;)) = (dz/da;)"i^/(z(x)), gives 

z^A'higz^) ^ ^ [j^a^" ~ i9^')') ■ (3-19) 



4 Random-matrix partition function for spanning forests 

Conside the set of pairs (G, F) where G is a connected graph with all vertices of degree k = h+2, 
and F £ J^{G). With our choice for the measure, the partition function Z{t,g,N) is given by 

(|2.9p . Call {Ta}a=i K{F) the components of the forest. Consider the edges E{F) C E{G) 

as marked, and the remaining edges as unmarked. Conversely, consider all vertices as marked. 
Now the quantity K{F) coincides with the number of connected marked components (possibly 
consisting of isolated vertices) . We say that an unmarked edge is an arc if it connects points 
on the same component, and a bridge if it connects points on distinct components. 

For each edge (ij) S E(G), call leg-decoration the introduction (in scries) of two interme- 
diate vertices, ij and j^: 

i j ^ i ij ji i 

• • • • • • 
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Figure 2: On top, a portion of a typical configuration of spanning forests on a 
random 3-graph, and a planar cubic tree corresponding to the right-most connected 
component. Below, on the left, a manipulation of the drawing which highlights 
the shape of the "effective vertices"; on the right, the diagram of the one-matrix 
partition function which contains the contribution of the original configuration, 
obtained shrinking the effective vertices to single points, by keeping memory only 
of the external-leg sequence. 



Call leg a decorated edge of the form {iij). This decoration provides a language for the 
mechanism of Wick contractions underlying the Random Matrix technique; the vertices coming 
from the expansion of the action determine the set of legs (with their cyclic ordering), while 
the choice of Wick contractions determines how the legs are connected. 

Leg-decorate all unmarked edges. For each component T of F, we define its border as the 
set of unmarked legs incident on T. As G is actually a 2-dimensional cell complex, and any 
tree is homotopic to a point, the legs on the border of a tree T have an induced cyclic ordering. 

Consider a component T with n vertices, together with its hn -\- 2 border legs, as a new 
tree T' with n vertices all of degree k, and hn + 2 leaves. We know from problem 3 that the 
number of such configurations is exactly Ah^ni and, if divided by the factor hn + 2 given by 
the cyclic symmetry, is A/, „ • 

Imagine to contract these trees {T^} to single vertices, without altering the cyclic ordering 
of the external legs. Only edges of the form [ij ji) survive, and the combination of all possible 
n- vertex trees, with n > 1, leads to an "effective" coupling c/„ for the resulting diagram, which 
coincides with G/F. Recalling that we have a factor t per component, and that we introduced 
the proper symmetry factor in the definition of we should write 

= tg-A' . (4.1) 



hn + 

The combinatorics of the resulting (fat)graphs G/F is then suitable for resummation through 
Random Matrix technique. Introduce a Hermitian matrix of fields M, and consider the action 



n>l 

= -\ix{M\l-2tA'^{gM'')) , 



(4.2) 
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this leads to our desired generating function: 

Z(i,.g,7V) = -^ln / dMe^^(*^). (4.3) 

A graphical explanation of the whole procedure is presented in figure [2] 

So, we still deal with a one-matrix theory, as in the case of pure-gravity, although, un- 
fortunately the "potential" is not polynomial in matrix fields (and this causes problems, for 
example, in the solution methoud through the resolvent function). Again we can reduce to 
eigenvalues, and obtain 

Z{t,g,N) = ^In |^dAA2(A)e-^^.^(^'); (4.4) 

V{X)^^a~2tA',{g\'^)). (4.5) 

where integration over values of A is intendend inside the analiticity region for 1^(A), that is, 
using the result of (|3.12p . 

\X.\<2^if^y. (4.6) 

We can reduce the potential to a polynomial, via the proper change of variable, inspired by 
relation ((37T9)) 

Xi{x,) = x,il - gx'l) ; (4.7) 
such that the interesting quantities change into 



d\,^dx,{l-g{h+l)x^); (4.8) 
|Aj - AjI = \{x^ - xj) - .g(a;f+^ - x'-^'^)] 

- \x^ -xj\-\l- gix'l + x'r^xj +■■■ + 4)1 



. , Jr-i^ , , ./n. (4.9) 



Vix) = y (l - 2.9 (l + TTTs) + (1 + t) x''^) . (4.10) 

Remark how the corrections to the measure and to the Vandermonde factor combine, to give 
the fcLCtor 

A,(f) = - 9{^^ + ^t'^s +■■■ + 4)) ; (4-11) 
such that now the partition function reads 

Z= ^ln^dxA2(f)Ag(x) exp(~iV^y(a:,)) ■ (4.12) 

5 A remark on the 0{n) loop- gas model 

Given the spherical 0{n) probability measure 

dfiia) = -^5{a^~l)d"a, (5.1) 

with n„ = , the lowest moments are given by 



(a") = (a W) - ; (aV) = -Sab ■ (5.2) 

n 
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Consider a graph G in which each vertex has degree at most 3: the action 

Sg^ H^ + nl3a,-S.j) (5.3) 

{i,])£E(G) 

is such that the polynomial expansion of the integrand in the partition function involves only 
the moments in (|5.2p . More specifically, the terms of the expansion of 

n ll + n[3d.,-dj) (5.4) 

(ij)6S(G) 

are in correspondence with the configurations a € ({0} U {1, . . . , where an unmarked 

edge {ij) (that is, = 0) corresponds to a choice of the summand 1 in the expansion of factor 
(ij), and an edge marked with colour a (that is, = a) corresponds to the choice of summand 

Integration over the spherical measure at each vertex leaves only with configurations of 
marked self-avoiding loops, weighted with a factor [3 per marked edge. Summing over loop 
colourings also produces the "topological" factor n per loop [10]. 

The same result would have been obtained for any spherical measure such that equations 
(j5.2p hold, i.e., up to a rescaling, for any spherical measure. We remark however that, if we 
choose to integrate the variables with the very special function such that not only equations 
()5.2p hold, but also all higher momenta vanish, the combinatorial 0(n)-model-loop-gas cor- 
respondence extends to generic graphs and actions, up to the possible appearance of dimers. 
Indeed, assume the action has a series expansion 

^ = EE/3fc"'('^^-'^j)'- (5-5) 

(y> fc>i 

As momenta higher than ((cr^)^) vanish, all the coefficients (3k with fc > 3 are irrelevant. We 
can equivalently parametrize /3i^2 as 

/3i=/3; /32 = -^/3' + ^7; (5.6) 

and the integrand of the partition function is 

e^H =\{(l + npY^ay^ + n^Y^^ay^A + R{a) , (5.7) 

(jj) \ a a / 

where the remainder term R{'j) is a polynomial in the fields, in which each monomial has at 
least either a factor af or a factor cr^(j\ with a ^h, and thus vanishes after integration, with our 
choice for the invariant measure. So we can consider a combinatorics of unmarked (summand 
1), marked with colour a (summand n^afa^) and doubly-marked with colour a (summand 
n^{af(jj)'^) edges. Variable integration produces 1 if all adjacent edges are unmarked, 1/n 
if two adjacent edges are marked, and with the same colour, or one edge is doubly-marked, 
and otherwise. So we are left with configurations of coloured self-avoiding loops and dimers, 
edges in the loops being weighted with a factor (3, and dimers with a factor j/n, Summing over 
possible colourings reproduces the "topological" factor n per loop, and rescales the weight of 
marked (but uncoloured) dimers to 7. The case 7 = in the action (|5.5p . on a cubic lattice, 
corresponds to the loop-gas problem studied in the literature, in [10] and subsequent works. 

6 Connection between the spanning- forest and the 0(n) 
loop-gas model 

The problem of counting configurations of self-avoiding closed loops, with a weight (3^n^ given 
to a configuration with £ loops of total length L, is expected to be a combinatorial variant of 
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the 0{n) model, in particular in the case fc = 3, where the combinatorial derivation is more 
transparent, and for this reason has been widely studied first by Kostov, and afterwords by 
many others, with a large number of interesting results [32-34,37]. 

Here we briefly sketck the derivation, in order to highlight the similarities with our partition 
function, in equation (j4.12p . It turns out that a stronger analogy emerges in the generalization 
of the loop-gas problem considered in section [51 so we will study the problem of counting 
configurations of self-avoiding closed loops and dimers, on random graphs of coordination 
k = h + 2, with a weight /3^n^7'* given to a configuration with £ loops, of total length L, and 
d dimers. 

One can study the problem with a random-matrix technique, introducing n + 2 Hermitian 
matrix fields, a matrix M for unmarked edges, a matrix A for the dimers, and n auxiliary 
matrices {i?a}Q=i, which mark the edges of the loops in one of n colors, in order to 
reproduce the "topological" factor n^. Thus we have the random- matrix partition function 

^i-d gas - ]^ In / d(Af, A, {i?4)e^'5(M,A,{iJ„}) . ^g^^) 



S = ti 



^(^M^ + A^+Y,El)+jM'' 

a 



(6.2) 



Oijl' 

If we perform the Gaussian integration of matrix A, we have 

S = tr( - \{M^ + y^El) + ^M^ + 7g''M2'=-2 ^ ^m""' E^M''~^' E^) . (6.3) 

a a,h' 

Consider the change of variables 

M UMJ-^ ; Eo, UEo,U~^ ; (6.4) 

where the Jacobian is just the one-matrix Vandermonde determinant, as the transformation is 
unitary on the auxiliary matrices Ea- Now, for any pair i < j, each term {Ea)ij appears in 
the action with a quadratic contribution proportional to 

{E^h{E^h{l - Pg'iXl + A^^A, + ■ • ■ + A^')) (6.5) 
and thus the Gaussian integration of E^ degrees of freedom gives a factor A-g/^{X) 

^i-d gas = ^^^J dxA'{x)A-;/\x) exp ( - ^ V{x,)) ; (6.6) 

y(^)-\^'-J^'-19''x^'-'. (6.7) 

If analitic continuation in n can be performed, at least for negative integers values of n 0, in 
the case n = — 2 we recover the partition function (|4.12p up to constants, with the parameter 
correspondence 

{ g'/g^l/l3^h + 2 + t- 

1 7 = -i(.g7.9r^(i + t)5-; ^ ^ ^ 

the only last delicate point being the fact that in section 3] we required all Xi to be in the region 
of analicity of the function A^ , while in this case we require that all the factors [1 — Pg' {h+l)x'^) 
are strictly positive, in order to make positive-definite the Gaussian integration of the auxiliary 



^Cfr. [36] for a discussion on this point. 
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degrees of freedom (clearly, l — a{h+l)x'l > for all i implies 1 — a X]^,' ^'i > for all pairs 

One should check that these requirements coincide. A one-line argument is that the 
quantity In A(2;) acts as a sort of "repulsive potential" from the border of the allowed domain, 
thus, as it coincides in the two cases, also the borders of the two domains must coincide. 

For h — I, the case of Nienhuis loop-gas, with no dimers, is obtained for 7 = 0, that is at 
t = —1. More generally, at arbitrary degree, a loop-gas model with vertex-disjoint loops and 
no dimers is still obtained for t = —I. We remark that this value of the coupling is special in 
the theory: the generating function of spanning forests is 'probabilistic' only for real positive 
values of t, by which we mean that each configuration takes a real positive weight, which can 
be interpreted as an (unnormalized) Gibbs measure. For t real negative, this picture does not 
hold anymore, in the combinatorial formulation in terms of forests. It does hold, however, for 
an alternate description, in terms of spanning trees only, weighted with some non-local factors 
related to "activities" , this being the original Tutte description of the generating function. 
In particular, for the generic Random Cluster model we have the definition of internal and 
external activities, while specialization to spanning forests gives unitary weight to externally- 
active edges (so that it is not necessary to count them), and gives a factor (l + t) per internally- 
active edge. So, in the interval t € [— l,-|-oo), the description of the generating function in 
terms of trees and activities is probabilistic. Tutte's notion of activity requires a choice of 
a linear order on the edge set (though the generating function of the activities is, in fact, 
independent of this order). 

Various other definitions of notions of activities exist, leading to the same Tutte generating 
function, in a non-obvious way. In particular, an astonishingly different notion has been recently 
introduced by Olivier Bernardi [38]. His notion only requires a cyclic ordering of the edges 
around each vertex (and, again, the resulting generating function is, in fact, independent of 
this choice). Not only it is interesting that, for a generic graph, this characterization strongly 
restricts the range of arbitrariness in the accessory ordering structure, but also, at our purposes 
here, that it is specially natural for graphs which are already embedded on the Riemann sphere, 
as is for the ensemble of "random planar graphs" 2-dimensional cell complexes arising from 
Random Matrix theory. In this case, the embedding naturally defines such a cyclic ordering, 
and no arbitrariness whatsoever is required in the Bernardi construction of Tutte polynomial 
(except for the choice of a single starting directed edge). So, this is a natural candidate for 
the construction of a probabilistic combinatorial expansion of our generating function, and 
the investigation of the model at t \ — 1, when internal activities are forced to vanish. This 
combinatorial approach could shed a light on the natural conjecture that the model is critical 
at t = — 1 for any graph degree, and in the universality class of Nienhuis loop-gas model. 

7 Perturbative expansion: spanning trees 

In the previous sections we sketched how to deal with the model at generic t with standard 
random-matrix techniques, even at finite N (and, for example, have access to higher-genus 
generating functions via loop equations). 

Beside this, a purely combinatorial approach based on the results of section |3] is enough to 
determine the planar partition function perturbatively in t at any given order, with a relatively 
small effort (also higher-genus quantities could be calculated with this technique, but this is 
not discussed here). This approach is interesting, not only by itself, and for the emerging 
combinatorics, but also for a comparison with the analogue perturbative calculations on flat 
2-dimensional lattices. 

In this section we perform the first-order calculation, concerning spanning trees, while in 
sections from |8] to [12] we will describe the higher-order technique. 

So, we deal with the problem of counting the pairs (G, T), with G a connected planar graph 
with V vertices, all of degree k = h + 2, and T a spanning tree on G, denoted by the symbol 
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T ^G.We have 



The graph G /T contams only one vertex, of coordination hV + 2. The edges form a hnk 
pattern connecting these terminations. 

The combinatorics of the coefficients, together with the constraint that hV is even (the 
ensemble of coordination-A; graphs with V vertices is empty if both V and k are odd, as there 
are no pairings of an odd number of legs) , produces 



Ai.yC{hv+2)/2 h even; 



Zi{g) = { (7.2) 
Z2 9^^'h,vC{hv+2)/2 h odd. 

V even 

It is worth stressing why A'j^ y is the appropriate coefficient: indeed it forces the graph to have 
at least one vertex (as q = our definition), and has an appropriate symmetry factor 
which accounts for the relative cyclic rotations of the link pattern and the leaves in the tree. 
So, in the case of h odd, scaling the index V in the sum by a factor 2, we have 



^^^y {2V{h+l))\ 

from which we have the asymptotics 



allowing to obtain the critical value of the coupling 

^eW = 2-()^^> (7.5) 

and the universal exponent —4, which is in agreement with the KPZ prediction. Indeed, we 
have just found that the string susceptibility 7 of the random-graph spanning-tree model, 
defined through the relation 

Z-^a;"n-3+^, (7.6) 

n 

must be 7 = — 1. On the other side, the central charge of the model in flat 2-dimensional space 
is c = —2. Finally, the KPZ relation consistently predicts 



^ ^ c-l-V(25-c)(l-c) ^^^^^ 

Furthermore, if we recall that spanning forests emerge as a limit g ^ of the Potts model, we 
also find that the value of the string susceptibility is in agreement with the formula of Eynard 
and Bonnet [39], valid for the Potts Model in the range q € [0,4], that, in the parametrization 
g = 2 — 2 cos(7ri') reads 

7-1-^. (7.8) 

From the ratio of two consecutive summands in (|7.3p . we deduce that the expression for Zi(g) 
is a generalized hypergeometric function of variable (g/gc)'^ . In particular, for the case of cubic 
planar lattices (shifting summation index from to — 1) we have the function 

^i(.9) = ^(2Fi(-|,-i;2;2V)-l) • (7.9) 
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In a similar way we can handle the case of h even. We have the formula 



V 

from which we have the asymptotics 

V 

which again gives the universal exponent —4, and the equation (|7.5p for the critical value of the 
coupling. From the ratio of two consecutive summands in (j7.10p we deduce that we deal with a 
generalized hypergeometric fimction of variable g/gc- In particular, for the case of coordination 
4 (also in this case shifting summation index from V to V — 1) we have the function 

^i(g) = ^(2Fi(-f,-i;2;33g)-l) . (7.12) 

As we are interested to the limit of large volume, it is instructive to analyze the approximate 
expressions (|7.4p and (17. lip . Parametrizing g/gc — e~'^ if h is odd, and g/gc — if h is 
even, and determining through Stirling expansion the overall constants in ()7.4p and (|7.1ip . we 
get the unique expression 

f U AA 

J- odd, 

fZi(g)^t^(/i)^e-2-n-4; ^^{h) ^ I ^jr^ (7.13) 

">i 2^1^^ /ieven. 

(The even/odd feature comes from the fact that, in the even case, the summands have a double 
density). The sum gives a polylogarithmic function, whose series expansion in e is 

Li4(e-2') = C(4)-2eC(3) + --- + ^e3lne+..- , (7.14) 

of which we highlighted the only term in the series which is not a pure monomial in e, but 
involves the non-analyticity Ine (this non-analiticity arising, of course, in addition to the fact 
that the series is not convergent if e < 0). We shall interpret the exact values of the terms in 
the series as being affected by the approximation of the summands for small n, however this 
general feature (that the first logarithmic term comes with a power e^), as it must be driven 
by the large-volume limit, is well captured by the approximation, and even the corresponding 
numerical coefficient 4/3 should be exact. Indeed, the (more painful) series expansion of the 
exact expressions (j7.9p and (|7.12p gives, for h = 1, 

1 /8192 16V27r\ / 14848 32\/27r\ 4„ , ^ 

and ioi h ~ 2, 

tK{h)~^Zi{g)= (^ip-24V37r^ + - 24^3^^ e + • • • + ^e^ In e + • • • . (7.16) 

So we conclude that the leading non-analytic contribution to Zi (in powers of e) is given by 



tZi[g) 



(0(e), analytic) + n{h)\e^\ne[l + 0{e)) 



(7.17) 
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8 Perturbative expansion: two components 



Our goal is to construct a perturbative expansion for the partition function (|2.9p . at every 
order in t, using the combinatorial results of section [31 The case of spanning trees described 
in section [7] was specially simple. Also the case of 2-forests is slightly special, for a number of 
concerns, while from the third order on we can give a general recipe (this is done in section [5]). 
Consider forests consisting of two trees, named as = (Ti,T2). We have 

'^'^ '^i-^f t,\p^^t{G/F)y ^^-'^ 

where summation over G is over connected planar graphs which contain the two trees Ti^2, 
and no other vertices beside the ones in the trees. The factor i is due to the symmetry 
under exchange of the two trees. Shrink the two trees Ti,2 to single vertices vi^2, and call 
deg(?;) the degree of vertex v. In order to have an integer number of vertices per tree, as 
-^1,2 = ■|;(deg(t;i,2) — 2), it must be deg(?;i^2) = 2 modulo h. 

Now we can deal directly with the graph G' = G/F, at the cost of a combinatorial factor 
for the number of "reconstructions" of the trees. The edges of G" are either bridges or arcfl 
For f > 1, consider the set 5*^ of graphs G' with i bridges. On any vertex, in the cyclic 
ordering, two subsequent bridge terminations are divided by a (possibly empty) link pattern 
configuration. The exact combinatorial factor is 

-jAh,NtAh^N2 , (8.2) 

plus the constraints that A^i and N2 are positive integers. The symmetry reasonings leading 
to this result are as follows. If we can fix a "canonical" leg among the ones incident on ui.2, 
we have a factor Ah^M, instead of the A'^ that one should combine with a summation on all 
the equivalent legs. A candidate canonical leg is a leg of a bridge edge. However, we have t of 
them, "equivalent" because of the cyclic symmetry of the definition, so we must sum over the 
marking of all the equivalent bridge edges, and divide accordingly by 

Thus, we have a generating function (in some parameters Zi\, for reasons explained in a 
moment) 



i\Zl Z2 



(8.3) 



and Z2 is recovered from this through the formal substitution z^ Afi.N, and z-^ ^ ^ 

if iV' ^ 2 modulo h. 

Remark the appearence of a recurrent combination 

q{z):^^^ = l{z-V7^4) (8.4) 

(this is the expression for the resolvent of the Gaussian theory), having the useful property 

q{z)^ = zq{z) - I . (8.5) 
Summing over i we are left with the function 

Q'{zuZ2) = J2 7(9(^1)9(^2))' - - ln(l - q{z,)q{z2)) . 



Recalling that the function A^^z) is such that 



(8.6) 



- z\A^{gz'^) - l)z-('"+2) = A,„(l (5„,o)ff" , (8.7) 



2Triz 



^This is said following the terminology of page [8] i.e. an edge is an arc if it connects points on the same 
component, and a bridge if it connects points on distinct components. 
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we can use contour integration to implement the formal substitution above. Define 

a{z) = z{Ah{gz'^) - l) , (8.8) 

then we have 



^ A — 1 O 



(8.9) 



We would now need to determine the leading non-analytic term in the limit g /* gc^ i.e. the 
two-component analogue of equation (j7.17p . For some reasons (namely, the cyclic symmetry of 
Feynman diagrams constituted of two vertices and £ edges in parallel, which leads to logarithmic 
series instead of geometric ones), this calculation turns out to be a subtle variant of our 
approach for forests with n > 3 components, analysed in the following sections, and is discussed 
in full detail only in section [TUl 

Here, and up to the end of the section, we will make a digression, in which we use a different 
method, w.r.t. the one of contour integration in complex plane, used in the rest of the paper. 
This tool succeedes in determining the nature of the leading non-analytic behaviour, but fails 
to provide the exact value of the corresponding numerical coefficient. 

We will not enter too much in the necessary mathematical background of the method (that 
we could call of "Levy calculus"), however such a pedagogical introduction should appear in a 
future companion paper, where a number of other statistical results are derived in this way. 

We go back to the understanding of Z2 as a sum over the number £ of bridges, and, given 
we have 2i link pattern configurations intertwined with the bridge legs, cyclically along the 
two vertices. So we have 



£,Afi,Af2>l 
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1^1 ,. . . ,z>2f ^0 a— 1 



(8.10) 



X ^ 



We can easily manipulate a number of overall factors. Define the rescaled Catalan numbers as 

Cn = 2-2"-iC„ , (8.11) 

such that the corresponding generating function has radius of convergence 1, and indeed these 
coefficients are normalized, X]n>o C*,! = 1. Define also the rescaled coefficients A^^jq as 

Kn - {2^^9c{h)YAn,N (8.12) 

so that also their generating function has radius of convergence 1 (it is however normalized to 
1 -f ■^). Through these coefficients, the leading exponential factors are clearly highlighted 

-, / XN1+N2 

Z2{g) = 2^ ^ -Au.N,A^.N2 - 

(8-13) 

X ^ ShNi+2,2(vi+---+v^)+t ^hN2+2,2{vt+i+---+V2i)+i W C'^c ■ 
1^1 ,. . . ,z>2«^0 a— 1 

The parameter i is the only variable which entangles the contributions coming from the two 
trees, so for a given value of I we can investigate separately each term 

^2(5) = E^/^(5)'; (8.14) 



i>i 



feig) ^ ^ ^ihN i—j ^hN+2.2{'^i+-+,.t)+e YlCi^^ ■ (8.15) 
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The quantity fi can be easily calculated exactly for £ = 1 and h = I, with a method analogous 
to the one of section [7l The result is 

fii9ce-n - E A.a.-iC.e-^'"-'^ = e/ (l - ^Fii-'j, i; 2; g-^^)) 



where we highlighted the leading contribution, and the leading non-analytic contribution. So, 
the leading contribution to fi is of order 1, while the leading non- analytic contribution comes 
from the cross product, and is of order Ine. If we factor out a term (i/e) in t'^Z2, we have 



a combination of the kind 

(^(e2 + ...) + B(e* + ---)lne) , (8.17) 

analogous in form to (|7.17[) . with the difference that the analytic term starts from a higher 
order. 

We will now study the function fi{g) for arbitrary values of and in particular in an approx- 
imation valid for large values. At this point, we make use of the core of Levy calculus, i.e. the 
generalized Central Limit Theorem for variables sampled from a heavy-tailed distribution (see 
for example [40, sec. 3.7]). Defining "the" Levy distributior|f| 

^cix) = .[^e-^x-l , (8.18) 

which is a special representative of the Levy family of alpha-stable distributions, for a = i 
and (3—1, with pseudo- variance c. We have 

e 

J2 Sn.u,+-+., n -C^{N + 0{e)) , (8.19) 

1^1 ,. . .,i^£>0 a — l 

and, in particular, typical values of N are of order l'^. The reason for this is that rescaled 
Catalan numbers have the asymptotics Cn ^ ^^^^ ^^^^ distribution on 

the integers C„ matches with Ci{x), and the sum of N independent variables with parameter 
a and pseudo- variances Ci follows an alpha-stable distribution with parameter a, and pseudo- 
variance c satisfying = ^^cf. A more subtle analysis (namely, a second order in Stirling 
approximation for the Catalan numbers) would allow to state in (j8.19p 

C^iN + Oie)) ^ Ce2(N ~ ^£ + Oil)) , (8.20) 

2 2 

SO, in order to control the errors deriving from the linear part in i, we will write N ~ b£ as 
argument of the Levy distribution. 

At this point an issue of factors 2 coming from an even/odd feature emerges. The point is 
the constraint that 2(z^i -t- • • • -I- I'e) + £ = 2 modulo h. If h is odd, this happens approximatively 
with a flat probability l//i, when N is large, regardless of £. If h is even, this can never happen 
if £ is odd, and happens approximatively with a flat probability 2/h, when N is large and £ is 
even. To be definite, we will consider the case of odd h, although it would not be difficult to 
treat also the other case. We get 

^ 1 

E w2,2(.,+...+..)+. n ^-^^ ^ (K^ - ^^)) ■ (8.21) 

i^i,...,z^£>0 a—l 



•^Some authors name this distribution after Levy, while others use the name of Levy distributions for the 
broader family of alpha-stable distributions. 
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Parametrizing g = gcC ^ , and using the asymptotic behaviour of A/i_jv (as deduced from (|3.9[) ). 
one has 

Af>l ^ ' 

We see how the integrand has a good scaUng for iV ~ £^ and e ~ with the correction term 
in h subleading (of order j). We wiU neglect it from now on. 

Approximating the sum with an integral (again, legitimate up to corrections of relative 
order j\ and using the known formula 

dx 



where Ki{x) is the Bessel K function of index 2, we get 



exp {-ax - ^) = ^K2{V^) , (8.23) 



Substituting into the expression (j8.14p for Z2, we have 



^ ' £>1 

In this case, the leading non-analitic behaviour in e must be deduced from the characteristics 
of Bessel function K2 near the origin (as the contribution at large I is suppressed both by the 
algebraic prefactor, and by the behaviour of K2 itself). We have 

i^2(2z)- (^^ + ---) - lnz(^y + •••) , (8.26) 

where the dots stand for further terms in a series expansion. So, the leading term in the square 
goes like while the leading non-analitic term comes from the cross product, and goes like 
—z^Xxiz K2{2z). Using also the formula 

^K2{x)^^^^\{l + 0{a'\na)), (8.27) 



we get 



and finally 



5^ ^K2{V~^er = (4C(7)e-^ + ...)-iln.(i + -- -), (8.28) 

e>i 



t'Z2{g) 



{0{e^), analytic) - ^^±^e'\ne{l ^ 0{e)) 



.29) 



The fact that the leading non-analytic behaviour is dominated by the cross product (thus 
breaking the symmetry between the two trees), and that the leading contribution to the sum 
(j8.28[) comes from small values of £, is a hint towards the fact that the generating function Z2, 
in its limit of large graphs, is dominated by forests in which one of the two trees is much larger 
than the other. We will come back to this point in sections fTOl and fTTl 

Bessel functions do not appear anymore in the paper, so there should not be confusion with 
an unrelated quantity, defined in the following, for which we use the letter K. 
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9 Perturbative expansion: higher orders 

Here we show how to construct a perturbative expansion for the partition function (|2.9p . at 
every order in t, using the combinatorial results of section [31 The cases of spanning trees and 
spanning forests with two components, described respectively in sections [7] and [51 were special 
under certain aspects, while from the third order on we can give a general recipe, being a 
slight modification of the technique of Cauchy integral explained at the beginning of section [51 
(conversely, we do not use anymore the tools of "Levy calculus" adopted in the remaining part 
of that section). 

At the light of the 'hardness' of the derivation for 2-forests, in comparison with the one for 
trees, the topic of the present section could seem an ambitious task. However, a more careful 
comparison of formulas (|7.17p and (|8.29p leads to more optimistic expectations. Indeed, from 
these partial results, it looks like the series Z{g) = ^^f^Znig), for g = gce~'^ near to the 
critical value, can be written as a double series in the variables t/e and e. In formulas (|7.17p 
and (|8.29p we highlighted both the tout-court leading behaviour, and the leading behaviour 
among the terms containing some non-analiticity for e — > (namely, a factor Ine). Of course, 
the more the latter is subleading, the more a refined control on the result is necessary, in order 
to extract a sufficient number of terms in a full expansion. Conversely, if the latter is leading, 
a first-order perturbative analysis should be sufficient. From the analysis of the cases with n 
components, n being 1 or 2, it can be conjectured that the disturbing analytic series starts with 
e", while the first non-analytic term occurs at order e'^lne (in agreement with the value of the 
string susceptibility, and the fact that \V\ ~ 1/e)- As a result, if such a conjecture did hold for 
all values of n, for n > 4 we would be in the easier situation in which first-order perturbative 
analysis suffices. 

So we start the analysis of the combinatorics for forests F = (Ti, . . . , T„), with n > 3, being 
spanning on a graph G y F. Shrink the trees to vertices vi, . . . , u„. The graph obtained so 
far has been called G' = G/F. Imagine G' as drawn on the Riemann sphere (i.e., on the 
plane, plus a point at infinity). Given an arc, two regions on the Riemann sphere are naturally 
identified. Say that the arc is contractible if at least one of the two regions contains only arcs. 
For the remaining edges, given two edges connecting the same pair of vertices (they could 
be two bridges or two non-contractible arcs), again two regions on the Riemann sphere are 
naturally identified. Say that the edges are multiple if one of the two regions contains only 
arcs CaU G"(G') the graph in which contractible arcs are removed, and G"'{G') the one 
in which both contractible arcs are removed, and multiple edges are replaced by single edges, 
by shrinking out the regions containing only arcs. This procedure is depicted in figure [H By 
construction, G'" does not contain neither contractible arcs, nor consecutive multiple edges. It 
is easily seen that, at fixed n, the number of planar graphs with these characteristics is finite, 
and the number of edges is at most 3n — 60 In a sense, the graph G'" describes the relevant part 
of the adjacence structure among components of the forest. We have a 'bridge' edge between 
two vertices if the corresponding trees are adjacent in the full graph, and a non-contractible 
loop incident on a vertex if the corresponding tree T has edges in G \ T incident on vertices of 
T with both terminations, and which leave some other components on each of the two sides. 

The removal described above corresponds to a resummation analogous to the one done for 
2-forests, in determining the combination l|8.6p . Consider the generating function for graphs 
G' with n vertices 

"■'-Sm n («) 

G' ' ^ iev(G') 

^Remark that, if n > 3, at most one of the two regions has this property. This is one of the reasons why the 
case n = 2 is special. 

^This is a consequence of Euler formula 1 12. Il l for connected planar graphs, i.e. with t) = and K = 1, 
V + F — E = 2. Furthermore, the absence of contractible arcs and of multiple edges implies that all the faces 
have at least 3 sides, so that 2E > 3-F, equality holding for triangulations. Solving w.r.t. F gives the desired 
relation 

E<3y-6. (9.1) 
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which we will put in relation with the desired perturbative partition function 

^" = E^^ (9.3) 

G' ' ieV{G') 

In equation (j9.2p . the contribution of an edge between vertices i and j is XiXj, so we can 
equivalently write 

<=E^^ n (-^-.). (9.4) 

G' ' ^ ^' fo-)e£;(G') 

Given a vertex with a certain number of terminations of bridges and non-contractible arcs, 
each interval between these terminations can be occupied by contractible arcs in an arbitrary 
link pattern configuration. Thus, we can "dress" the edge terminations with the substitution 

Xi^XiC{xi), (9.5) 

and then restrict the sum to all graphs G" which do not contain any contractible arc, 



1 " 



(9.6) 



e-E(G') 



Then, multiple edges can be resummed. For each €-uple of multiple edges between vertices i 
and j, we have a contribution {xiC{xf)xjC{x^)Y . As each edge in G'" can be originated by 
shrinking an arbitrary number £ > 1 of multiple edges, we can replace them with dressed single 
edges through the substitution 

x,U[xJx,C(x^) ^ ^_^^c{x^)x,C{x]) ' ^^-^^ 

to be performed in the second formulation of quantity (|9.6p . together with a further restriction 
in the sum to graphs G'" which do not have multiple edges 

7' = V ^ TT a:»G(x-)x-,G(a;j) 

" Z.|Aut(G"')| \i-l-x,Cixl)x,C{x])- 

£E{G"') 

We claim that the formal substitution rule x^^^'^ A^^jsig^ if TV is a positive integer, and 

— > otherwise, applied to Z^, gives Z„. The only delicate point is the symmetry factor 
pertinent to vertices, i.e. we need to prove that a canonical way exists for marking a single 
leg incident to any vertex, without any other symmetry factors. Label arbitrarily the vertices 
with integers 1, . . . , n. This accounts for the already included obvious factor |Aut(G"')|~^. 
Remark that, as G'" is connected, each vertex must have at least one incident bridge, and that 
a canonical spanning tree on G'" is easily constructed, e.g. the depth-first tree starting from 
the 'root' vertex 1. 

Note that this depth-first tree has a natural orientation towards the root, so that any non- 
root vertex has a single outgoing edge, oriented towards his neighbour with smaller index. 
Also, this tree has a natural planar embedding in H (with the root vertex on 5IHI), and thus a 
natural oriented loop which encircles it (see figure [3] for an example) . 

For any non-root vertex i, having vertex j as a neighbour with minimum index, we will 
choose as "canonical" leg one of the bridge edges between i and j. We have in general L > 1 
such edges. As n > 3, at least one among i and j has other neighbours, so that the i-uple 
of bridges is univocally splitted into m €-uples of multiple edges, and in any ^-uple a "first" 
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Figure 3: On the left, a typical diagram, with labeled vertices. Bridges are drawn 
in bold. On the right, the corresponding depth-first tree, on the graph in which 
non- consecutive multiple edges are replaced by a single edge, and arcs are dropped. 
Remark that, in building the depth-first tree, the planar embedding of the diagram 
is discarded, and a new one is induced by the vertex labeling. In particular, the 
loop surrounding the tree is the contour of the light-gray region. 



edge in a cyclic ordering is identified (this is a difference with the n = 2 case, where no "first 
leg" can be defined). We will use as a canonical leg one of these first legs in a -^-uple, so we 
are left with m > 1 choices. We have to prove that a canonical choice among these m can 
be performed. At this point, the loop surrounding the tree plays a role. Indeed, the chosen 
pair (ij) is traversed by the loop, from i to j, exactly once, and thus a single "next" vertex fc, 
such that (jk) follows {ij) along the loop, is identified. Because of planarity, this vertex must 
be contained in exactly one of the regions identified by the m multiple-edges, and thus it also 
identifies a canonical choice for a leg, as was to be proven. Now that we have many reference 
points, for the root vertex any reasonable choice makes the game. For example, one can use 
the other termination of the edge containing the canonical leg of the neighbour with smaller 
index. This completes the discussion of the involved symmetry factors. 

So, we can use the technique described in formula (|8.7p . in order to extract Z„ from Z^, 

i.e. 




(9.9) 



The quantity in (|9.7[) . when expressed in terms of inverse parameters {z^ ^}, corresponds to 



q{zi)q{z2) 



■Q (Z1,Z2), 



(9.10) 



1 - 



1 



q{zi)q{z2) 



Zl Z2 



where q{z) has been defined in (|8.4p . and Q 
For future convenience, we also define 



(^1,22) is a useful combination. 



(3+(zi,Z2) := l + Q {zi,Z2) 



1 



(9.11) 



1 - q{zi)q{z2) 



and remark that these quantitcs can be restated as 




(9.12) 



(9.13) 
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Figure 4: On the top, a typical configuration of graph G' with n = 4 vertices. 
Contractible arcs are drawn as thin edges, and dotted lines collect multiple edges. 
On the bottom, the resulting graph G"'{G'). 



and in particulajfl 

Q-(z„z,) = ^P=; (9.14) 
Q"(....) = -^^^^^^^. (9.15) 

Zj 

In equation (|9.9p . expanded with (|9.8p . we can interchange integrations and summation over 
the graphs. Then, each term of the sum can be interpreted as a (coordinate-space) "Feynman 
diagram" : indeed, we have the appropriate symmetry factor, one integration per vertex, with a 
proper "measure" factor a{zj), and a product of "propagators" , or Green functions, Q~ {zi, Zj) 
corresponding to the edges of the diagram. 

G'" 3 = 1 e={i,j) 

eE(G"') 

In section [TU] we wih deal with the problem of understanding the appropriate contours for such 
an integral, and we will prove that a valid contour exists if and only if <? < gc{h). The integrals 



^Equation l|9.12|l comes from Il9.13| l through a I'Hopital Umit, while Il9.13| l. which coincides with 119.15^ 
trivially, is related to the definition I I9.10I I through the use of property 118.51 1. Namely, calling q = q{z) and 
q' = q{z'), by verifying that 

qq' ^ q- q' 

1 — qq' z — z' 



23 




Figure 5: Diagrams of the perturbative expansion for three-component (top-left) 
and four-component forests. Linear combinations, leading to the appearance of 
Q"*" propagators, are chosen in order to minimize the total number of diagrams. 
The coefficients, deriving from this combination of symmetry factors, are indicated 
next to each diagrams. Remark that the first non-contractible arc appears at order 
3, in the right-most diagram, while the first non-consecutive multiple edges appear 
in some of the diagrams at order 4- 



for the third and fouth order are 

3 



Z3= ^ 



n 27ri g'9i2'923('93i + 



2m 



2Ql2'923'534 ( l^Qt^Qti ''^ Q22Qt'iQ2'i 



(9.17) 



(9.18) 



1 (+2 1+3++ -'-+ + 

~^ '^Qi2Qi'iQii\Qii '9i2 + gQii ^ QiiQi2Q2'i^ '^Q2iQn^'^Qi2^ Q2i. 



where short notation Q^- stands for Q^{xi,Xj), and in many cases the contributions of more 
diagrams are collected together, with the combination = 1 + Q^. The related diagrams 
are shown in figure [51 where dashed lines denote propagators Q J- . 

Here we remark that this approach in some sense "captures" the degrees of freedom of the 
perturbative theory: we sum over the contribution of an infinite number of degrees of freedom 
(as we are summing over graphs with an arbitrary number of vertices), and we end up with 
a finite diagrammatics, i.e. a finite sum of finite-dimensional integrals. Remark in particular 
that, as expected, the diagrams in our family are all planar. 

A detailed discussion on how to practically estimate the result of these integrations appears 
in section [TUl 

We remark that most of the techniques outlined in this paper for the spanning-forest model 
immediately generalize to variants of the model, in which more general weights are chosen for 
the trees 



iV-2t)(G) 



\V{G)\ 



|Aut(G)| 



F^G l^^^^'^/^)! T„eF 



(9.19) 



here the weights w{T) for the single components (i.e. the trees) of the forest only depend on 
the structure of the tree (and not on the embedding into G). Our case is u'(T) = 1. Other 
interesting cases are 'w{T) = \T\ {rooted forests, cfr. below); the generalization u'(T) = [Tl"; the 
case w{T) = X]i>e\/(T) '?dcgy(t>)i which is the generating function for the distribution of vertex 
coordinations. In all these cases, the weight w{T) reflects into a modification of the generating 
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function a(z) in (|8.8[) . while all the formulas in which a{z) appears implicitly are valid in 
the generic case. In particular, the case w{T) — \T\'^, for various values of i', explores all the 
asymptotic behaviours of a(z), "shifting" the spanning-tree string susceptibility 7 = —1 to 7+;^. 
This is absolutely trivial in the case of trees, as |r| = |F(G)|, but has substantial consequences 
on the measure for spanning forests at t > 0. In particular, for j/ = 1 we have a model of rooted 
spanning forests, a considerably easier variant of spanning forests, corresponding to a massive 
perturbation of the pure graph-Laplacian implicit in Kirchhoff Matrix- Tree theorem. 



10 Evaluation of diagrams 

In the previous section we stated some combinatorial quantities in terms of contour integrals, 
which formally express certain convolutions of generating functions. In this section we discuss 
the analytic prescriptions on how to perform these integrations, and give a general technique 
to reduce our integrals to integrals of real- valued functions over real intervals. 

We give a number of examples, in increasing order of difficulty, with the aim of gradually 
introducing a set of tools, more extensively used in the following sections. For some of these 
examples, also the results are of direct interest. 

The easiest example of this method is still given by Catalan numbers. Suppose one wants 
to deduce formula (|3.3p from the generating function (|3.2p . Then, one can perform a contour 
integration 

C„./Ji£-1^^,-^ ,10,1) 



2TTiz 



It is legitimate to deform the contour of integration, up to encircle the whole plane except for 
the cut [i, +00]: 




For n sufficiently large (n > — ^ suffices, so for any integer n), the integral on the large circle 
vanishes, and we are left with the integral on the two sides of the cut. For z = x±ie, with x real 
larger than i and e real positive infinitesimal, one has -^/l — 4:{x ± ie) — ^i\^/4x — 1| + 0{e), 
and one can write 

which gives indeed Catalan numbers. 

A next case, of intermediate difficulty between the previous example and the general- 
diagram case, is the case of spanning trees on cubic lattices. Suppose we want to deduce 
the result of equation ()7.9|) with no use of the explicit coefficients C„ and A'^ „, but only 
through the generating functions C{x) and A[{x). Start consider the case of generic h. Call 
n the index for the tree size, and m the one for the number of edges, then the spanning-tree 
partition function reads 

Zl (g) = 9"K,nCmS2m,hn+2 ■ (10.3) 
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In generating function, 



„hn+2 



C{a 



-2m 



and contour integration can be used in order to reproduce the delta constraint 

dz 



Zi{g) = j> 



2-Kiz 



For the case h — \, using equations (|3.2p and (|3.17p . we have 

-1 + 6gz - ^{gzf + (1 - ^gzY- 



z'A'.igz) 



-C{z-'^) = 
z 



z — vz^ 



(10.4) 
(10.5) 

(10.6) 

(10.7) 
(10.8) 



Now it is clear which prescription we should use for the contour integral: we should encircle the 
origin, being in the radius of convergence for both series C{x) and A\{x), thus 2 < |z| < \/(Ag). 
As expected, no contour can be found when g > gc = 1/^. In terms of Cauchy integration, 
this prescription states that the contour should encircle the cut going from —2 to 2, and leave 
outside the cut going from \/ [Ag) to infinity. If we deform the path in order to have contribution 
only from the [—2,2] cut discontinuity, as in 




we have ^ 

^'^■"^ = / 2 £ y4^( -l + dgx- &{gxf + (1 - Agx)^) . (10.9) 

It can be seen that this is the appropriate expression. Indeed, if we expand the generating func- 
tion A'i{x), keep only even orders (because of the parity of the other factor in the integrand), 
and use the formula 



f ^^^^^fl—i^Cn (10.10) 



we obtain the series coefficients ()7.3p . 

What would have happened if we wanted to evaluate the integral by the change of variables 
which makes J^](x) an analytic function? The new variable would have been x' = a;(l — gx). 
For g infinitesimal, the cut [—2, 2] moves infinitesimally to the two new values 



l-^^r^^2 + 45+--., (10.11) 
~ ^-^/^ = -2 + 4.g + ..., (10.12) 



2ff 
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Figure 6: The counter-images of the interval [—2, 2] under the map z = x{l—gcX ), 
in the complex plane for x. Here h = 7. The bullets correspond to the counter- 
images of the origin. The dashed circle denotes the largest disk centered in the 
origin which does not intersect any of the external cuts, and is thus the radius of 
convergence for the integrand corresponding to a typical diagram in the Feynman 
expansion. A valid contour must encircle the internal cut, staying within this disk. 



and a new cut, to be left outside the contour of integration, appears, between the values 

2-4g + ---, (10.13) 



l + x/T~~% 1 



25 g 

.ii-^ = i + 2-4, + .... (10.14) 
2.9 9 

Above the critical value gc = 1/8, the two solutions and x'j^, instead of being radially 
ordered, are complex conjugate, and there is no contour in which the convergence of the 
generating function is assured. 

This phenomenology occurs also at higher values of /i, with a slight even/odd difference. 
Indeed, one can study the counter-images of the cut z S [—2, 2] under the map z = x{l — gx^). 
The value z = has /i + 1 counter-images x = {0} U {g~'s^e^^}o<£<;i, which 'mark' the h-^ 1 
counter- images of the cut if g < (where the critical value gc is defined in (|7.5p ). However, 
at g = gc, two real counter-images of z = 2, for the central cut and the cut with label £ — 0, 
collapse, at the point x^, = 2-^iii. If h is even, by symmetry this also happens to two real 
counter-images of z = —2, the ones for the central cut and for £ = h/2. An example for h = 7, 
at the critical coupling, is illustrated in figure [6] 

Another check of the validity of equation (|10.9|) . and a different approach to the result of 
equation (|7.17p . comes from deducing that, for /i = 1, 

Ziig) = (0(1), analytic) + K(l)^e3lne(l + 0(e)) . 



This is equivalent to state that 

Q3 



^^3^1(9) = 8K(l)lne + 0(l). (10.15) 
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Indeed, in (|10.9[) the parameter e occurs only in the "measure" A[{gz), and one has, for z ^ 2 
and g = gc<i^\ 

TT-^iz (qz)) ^ — 4(7z 2 -|- more ree. terms 

V )) 8 1252 ^ y ! s ^^^^^^^ 

= —2(1 — 4gz)^2 ^ more reg. terms. 



The semicircle factor ■\/4 — z^, near the singularity z = 2, is approximated by 2V2 — z. Then, 
integrating in some window [2 — a, 2] near the singularity, we get 

Zir^^- [ dzV2^fl--e-'')~' lne + C'(l) = 8K(l)lne + C'(l), (10.17) 

J2-a V 2 / TT 

as was to be shown. 

Now we repeat the procedure above in the case of ^2(3), in order to determine the coefficient 
of the leading non-analytic term, at least for h = 1, which was lacking in section [5] From the 
Levy calculus arguments of that section, we have an expansion of the form 

Z2{g) = (0(1), analytic) +ce2 In e(l + 0(e)) , 

with c an unknown numerical constant. This is equivalent to state that 

l^Z2ig) = cAne + 0{l). (10.18) 

We have to go back to the expression (|8.9p . This expression is symmetric under exchange 
zi ^ Z2, and the parameter e occurs only in the "one-body measures" a{zi) and 0(^2)- So we 
can write 

2de^ 2de^2J 2^ J 27ri ' ' 

1 f dzi / \ fdz2a{z2) , 



27ri \de^ ' 'J J 2m ' ' (10.19) 
=: + 

We will prove in the following that the second summand, i3, does not contribute to the leading 
non-analytic part of Z2- We start by calculating the leading expression corresponding to the 
first integral, A. It is useful to integrate by parts w.r.t. variable zi, in order to get a restatement 
which avoids logarithmic functions. This is affordable in principle also for arbitrary h, because 
the function a(z), of which little is known, has however a simple primitive (cfr. equation (|3.18p ). 
and would give 

^ ^ / 1^ (£(^MU..l'))) ^ / M . (10.0) 

The two derivatives acting on variable zi produce a singularity ^ 1/ z near zi = 2, up to more 
regular terms, so that, neglecting other more regular terms, we can replace zi — s- 2 in all non- 
singular expressions, and in particular in Q", thus factorizing the two integrals. As discussed 
later (in equation (jl2.6p ). the remaining integral in Z2 is regular for e ^ 0, and named Kfi. So 
we can write 

Kh f Az f \ 1 

-4 = ^ ^ — (^^(z2AU5^''))j ^7p=^ + i^O'^e reg. terms. (10.21) 
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so that 



Again we can deform the contour in order to take contribution from the cut, and get 

-4 - 1^ y ^ dz (^^{z'A'^igz'^))^ -^=L= + more reg. terms. (10.22) 

and one has, for /i = 1, z ~ 2 and g ~ gc&~'^ i 

— {z'A^{gz))^--^{l-Agz) ^ + more reg. terms ^^^^^^ 
= 4(1 — 4gz) ^ + more reg. terms, 

C dz^^ fl-|e-^r^ = -— Xilne + 0(1). (10.24) 

J2-a V2 - Z V 2 / TT 

More generaUy, for arbitrary h, going back to (|10.2ip . we have 

^iz'A',{gz')) = ^ [z^^^ ~ 3z^ + 4) {7?^^{gz^)) . (10.25) 

Higher derivatives lead to stronger singularities, so the leading order is due to the second-order 
derivative, and we can write 

^{z^A'.igz'^)) ^ l^z'^iz'A',,igz'^)) = ^z^^-^z) . (10.26) 

Having replaced derivatives w.r.t. e with derivatives w.r.t. z allows us to perform the g- 
dependent change of variables which simplifies the measure. In changing variables from z 
to a;, the Jacobian in the remaining derivative and in the measure simplify, and we have in 
general 

it n.) = / i| ( f (.(1 + . (10.27) 

while in our specific case 

^J {x{l + gx'^)) -4 

We can now reduce to an integral over the cut discontinuity due to the square root at de- 
nominator. Suppose that h is odd, so that we have a single point of singularity (instead of 
two symmetric ones). This singularity occurs at the rightmost extremum of the central cut 
(cfr. figureini), where two counter-images of z = 2 coincide for g /" g^. Denote with index 

1 the most internal solution, and with = 2-^^^ the limit for g = gc- 

As we are only concerned with the leading singularity, which is due to the square root at 

denominator ^ (a;(l -f- gx'^)) — 2, in all other factors we can replace the value of g with g^ and 
of X with x^, up to less singular terms, and write 

A - ^ r d. . 1 . _ , ^'^^ . = S r . . / . • (10.29) 



Vl^(x{l + gxh))^2 ^^"^ J sjx{l + gxf' 



Furthermore, in the polynomial of order h+1, x{l+gx^) — 2, we should highlight the two main 
roots X = xX' , and for the rest we can replace the other roots with their values for g = gc- 
Using the fact 



x{l + gx^)-2 _ Id' 



i2 



lim ^ ^= lim - — (x{l+gx'') 



(10.30) 
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we can write 



^ , ^ ^^^( - .n^- + Om) , (10.1, 



tt/i^ Jo h y/x{Sx^ - x) 7r/i3 

with the shortcut Sx^, :— x?-* — xi^K As we have 



we have 



so that we get 




(10.32) 



(10.33) 



■4^- ^^\l+^ (hie + 0(l)). (10.34) 

This is in agreement with the special case ft, = 1 of equation (|10.24p . 

For what concerns the integral B, we have that the two one-body measures produce a 
singularity ~ fo'^ of the variables 2:1^2 2. Beside this, there is only a two-body 

interaction factor Q'{zi,Z2), which is regular unless both varables approach 2, in which case, 
for zi,Z2 along the cut [—2, 2] and near 2, it behaves as 

Q'(zi,Z2) - iln(max(2-zi,2-Z2)) (10.35) 

so that, even taking directly e = 0, we get a finite result after integrating near to the only 
potential singularity 

Boc [ dzi [ dz2 = 2a(lna- 1). (10.36) 

Jo Jzi 

By collecting pim)l . PUM)) and pU^S]) we get 



(0(1) analytic, 0{e^)) , „ 



Vh+lKh 2, 
— ; e m e 



(10.37) 



The integrals arising from the full perturbative expansion are similar, although more involved, 
and, of course, multi-dimensional. However, as from three components on the symmetry factors 
are easier to handle, there is a small variety of fundamental ingredients: only the one-body 
measure a{z) and the two-body interaction Q~{z,z') occur. Consider a diagram D with n 
vertices and edge-set E{D). The integral is of the form 

^(^)=nf^^^ n (10-38) 

Consider the change of variables Zj — Xj{l — gx^). The measure changes as 

dza(z) = Ax 5a;''+^(l - {h + l)gx''). (10.39) 

We will keep using Zj as a shortcut to Xj(l — 52;^), and qj as a shortcut of q{zj), when this does 
not cause confusion. Recall that it is the expression z | — 4 in q{zj) which is discontinuous at 
the cut, and that, in the complex plane for x, the cut to be encircled by the integration contour 
is the real segment [x-{g),x+{g)], with < 2^^, and x+{g) for g — + gc- 

There are two 'easy' contour prescriptions, in the complex planes for Xj's, which are valid 
for any value of h> 1, and g < gdh). One is the circle of radius 2-^iii, so that one could use 
an "angular" parametrization Xj = 2-^iiie*^^, with 9j G [0, 27r]. A second one is to consider 
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a "cut" integral, i.e. integrate along the sides of a rectangle, centered at the origin, and of 
half-sides 2-^^ x 5, in a limit S ^ 0. In this limit only the cut discontinuity survives. 

In section [TT] we concentrate on the angular parametrization, when trying to extract the 
leading behaviour for g ^ gc- Here we discuss briefly the way in which the cut integral should 
be performed. Concentrate on a given variable Xi. Through the expressions (|9.14p and (|9.15p . 
the product of edge terms is, up to a prefactor not involving Zi, a function of the form 



n 



- qj 



with m's being non-negative integers. The polynomial in qi in the numerator can be reduced 
to a linear function in qi, of the form Po{zi) + Pi{zi)qi, by iterated use of (|8.5p . If rrii is even, 
only the square root cut discontinuity caused by the ^ zf ^ 4 term in qi gives contribution to 
the integral, so that we have 

dzj Po{zi) + Pi{zi)qi 



TTdz," ■ ■ ■ (p -—^a{z.,, 



^ (10-40) 



Azj a{zi) Pijzi) 



If instead rrii is odd, a cut term is already included in the denominator, so only the regular 
terms in the numerator contributes 



Poiz,) + Pi{zi)qi 



(10.41) 

dzi a{zi) 2Po{zi) + ZiPi{zi) 



2^ ^4^ (4 - zf)"^ UU^^ - ^.Y' 



In particular, a simple case of these integrals, in some variable Zi, is when i is adjacent to only 
one vertex j, through £ edges. In this case we have 



J 2^1 \ z — w 



i_2 27r ^ "\z-wY ^\2h+lJ V"^ '2) ^ ' 



(10.42) 



h 

The case ^ 1, in the specialization h = 1 and g = gc, gives 

^CiU \ I '^^ t xQ{z)-q{w) dz V4- 



(10.43) 



= — 2 - 3w) + 1 - u; - 4 gH o " 

37r TT + 2 V w + 2-2 

For w near to the endpoint of the cut, 

<S>^^\2 + x; gc) - f 3 - - 2^^- -xlna; + — (21n2 - 1) .t + 0(a;i) . (10.44) 

\ on J TT TT 

A calculation similar to the one performed in p0.43p leads to -^(I>(^)(w;gee"^)|^ _^, whi ch is a 
long expression that we do not write here. However, we report the equivalent of (|10.44p 



^$W(2 + a;;gee-^) 



e=0 



2^/2 ( A^/2 r- \ 1 

-^lna;+ 4+^f 6\/21n2 + ©(a;^ ) . (10.45) 

TT \ 6TT 
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Collecting the two results of (|10.44[) and (|10.45p we finally obtain 

$(i)(2 + a;;gee-^) = ( 3 - - 2^^ - -x Inx + — (2 In 2 - f ) x 

\ 37r / TT TT 

2/2 / 4/2 \ ^^"-^^^ 

+ — elna;+ 4+^— -6y21n2 e + 0{xh,x^. 

We thus see that the limit for w — *■ 2 of <i>'^-' is finite. Curiously, it coincides with the one for 
$(2), that is 

<i>«(2;.9) = <i>(2)(2;g), (10.47) 

as the difference in the integrand is given by a factor (cfr. the general expression (I10.42p ) 
(i) '^^^w-z^^ ' which goes to 1 identically in z, for w ^ 2. 

11 Non-analytic behaviour for g Qc 

In order to highlight the leading contribution to {Zn}n>3 in the large-volume limit g Z' gc-, 
we consider the series expansion described in section [HI in "angular parametrization" , that is 
parametrizing the 's and g as 

x,= \ ' e''^; 9 = 9ce-'. (11.1) 

In particular, the recurrent combination gx^'' becomes -f^r[^~'^'^^^^ ■ The various relevant quan- 
tities become 

= I e-+^('^+^)^(l - e-+-«) =: ^,{9) ; (11.2) 

z = 2e*^(l + i(l-e-^+'''"')); (11.3) 



- 4 = z^l - e-2'e(i + i(l - e-e+«/ie))-2 . (114) 

We will also adopt the shortcut f{9i, 9j) for Q^{z{9i), z{9j)). 

There are various potential sources of non-analiticities for g y gc-, due to non-regularities 
of the integrand for 6*^ ^ or 0^ — 9j — > for the various indices, when e approaches 0. 

A first singularity may come from the pole Zi — Zj in the denominators of Q^{zi, zj). 
However, for generic values of Zi , this is not a true singularity, as, for Zi — !■ Zj , also the numerator 
vanishes with the same behaviour (this is what allowed us to determine an expression for 
Q^{zi, Zi)). A first source of true non-analiticity comes, in the limit 9j — *■ O0from combinations 
of the form 

1 - e^'+''''' =e-ihe + \h^9^ + ©(e^, e9, 9^) . (11.5) 

Consistently with the fact that we drop higher orders in 6, we may adopt this approximation 
only in some small window 9 G [—J, (5], 5 <^ 1. 

We expanded at second order in 9 because, in the integration, only overall even monomials 
contribute. The measure ^(9) gives at leading orders, besides a factor as in (|11.5p . 

^_,+,(h+2)0 ^ 1 _^ ^. 2)61 + 0(e, 6*2) . (11.6) 
The product of (fTT3)) and pXB)) gives 

e-zh9+ ^^^^^^\ ^ + 0{e\e9,9^). (11.7) 



^Or 9j — » TT, if /i is even, however we neglect this for simplicity. It would be easy to reintroduce certain 
factors 2 overall at the end, in the case of even h, while the present treatment covers the case of h odd. 
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We will see in a moment that terms odd in 6 do not play any role in this measure. 

Indeed, a further source of non-analiticity for — > is the combination 2 — 2, appearing as 
a factor in the square roots inside terms . In this case, a stronger cancellation, also of the 
terms linear in 9, occurs. We get for the combination in (|11.3|) 

z = 2 + 1 (6+^^02^ +0(e2, 60, 03), (11.8) 

and thus we can rewrite the expression in (|11.4p as 

^^^ = l[e+'^^e^)+0{e\ee,e^). (11.9) 

This proves that the expressions Q^{zi, Zj), both in the case i = j and i j, are even in dj at 
leading orders, so we can drop out odd terms in the one-body measure d6jfi{6j), when \9j\ < S. 

This analysis holding in an interval of small 6 does not mean that the remaining part of 
the integral is negligible, and in general this is not the case. However, we can keep control on 
our expressions by considering an exact subdivision of the contribution of a diagram (jlO.BSp . 
Indeed, in general, for a contour integration on a path 7, and two points a and b on it, we can 
write f^dzf{z) = dz/(z) + J^^ij_^^-^ dzf{z). For our angular integrations, we can thus 

distinguish between the small-^ and large-^ behaviour by choosing to decompose the (periodic) 
interval [0, 2tt] into [—6, S] and [6, 2tt — 6], and denote by V C V{D) the set of angular variables 
for which we integrate in the first interval: 

I{D)^ I{D-V')- (11.10) 

V'<ZV(D) 

Remark that, while each T{D; V') in (jll.lip is a function of S, the sum T{D) in (jll.lOp is 
independent from its value. This arbitrariness will be exploited in the next paragraphs. 

Both the one-body function fJ,{Oj) and the Green function f{9i,9j) are of order 1 if the 
corresponding angles are larger than S. In particular, X(Z); 0) gives a contribution of order 1, 
analytic in a neighbourhood of e = 0. Thus it should be neglected even if it were the leading 
summand, in a way similar to what was shown to occur for the 1- and 2-component cases. 

The function n{9j) is given by (|11.7p for 9j small, and thus gives "small" factors, of which 
we can give a dimensional estimate 



)r^6ie + 6^). (11.12) 



Factors f{9i,9j) occur for non-contractible arcs, with equal indices, and for bridges, with 
distinct indices. In the first case, if 9 is small we get the expression in (|11.9p to the power — i 
(plus a subleading contribution of order 1): 

f{9„9,) ^ 4f ; l,^,, - + ^T'^ ■ (11-13) 

In the case of bridges, we have two non-trivial cases, depending if only one angle (say, 9i) is 
small, or both angles are small. In the first case, one can neglect the contribution of the small 
angle in the combination, as it does not produce any leading singularity, and write 



fi9.,,9,) ^ i [^|0±il ^ 1 ( jfl±|±i ) il + 0{V-e,S)) (11.14) 




where the leading expression is of order 1 by our previous assumption that j 1^ V . 
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The result above seems to suggest that, if instead both angles are small, a singular behaviour 
may arise, similarly to (|11.13[) . Indeed, in this case, we have for the leading part in f{9i,9j) 



Iv/if^-yipi 1 v/4(z, - 2) - ^A{zj - 2) 



2 Zi- Zj 2 (z, - 2) - [zj - 2) V^T^ + 

so that we obtain 



(11.15) 



^ ^(6 + ^^)-^. (11.16) 



/i V ^^ 2 '^i ^ Y ^ 2 '^i 

Remark how the expression (I11.13P is a special case of (|11.16p . with no need anymore of 
I'Hopital limit. At the end, we get 

X(Z?;r)~ (<5(e + 52))l'"l(e + ,52)-i.#{fe)GB(i3):..jev'}^ 

The expression above is just a crude "dimensional" analysis. However, it allows us to under- 
stand which terms, at fixed n and in the double sum over D and V ^ are dominant in the e — > 
limit. At this aim, choose -^/e ^ i5 ^ 1, and call E'{D, V') the set of edges in the subgraph of 
D induced by V . Then (|11.17p becomes 

We call order of a pair {D,V') the integer \E'\ — i\V'\. Pairs with maximum order give 
the leading contribution. After some reflection, one understand which families at fixed \V'\ 
minimize the order (besides the trivial case \V'\ =0, which gives order 0): 

Proposition 11.1 All and only the pairs {D, V) minimizing \E'\ — 3|t^'| at fixed n and h = 
\V'\ have the following defining properties: 

• The diagram D is a triangulation; 

• All vertices in V{D) \ V' have degree 1. 

Remark that the converse of the second claim, that all vertices of degree 1 are in V{D) \ V' , 
is not true. However, this case leads to terms which are strongly subleading in the fixed-n 
(and arbitrary h) analysis, and trivially reabsorbed in the contributions from smaller values of 
h. So, with an abuse of definition for the class of leading diagrams, which is justified at the 
light of the forthcoming equation ()12.3p , we will restrict our attention to diagrams in which all 
vertices of degree 1 are in V{D) \ V' . 

For these diagrams it follows easily that, for all the pairs {D,V') as above, 

• The subgraph D" C D induced by V\ and with loops dropped out, must consist of 
triangles and consecutive multiple edges; 

• For each pair of consecutive multiple edges in D" , incident on i and j € V , a non- 
zero number of vertices not in V' must be attached to i in the cyclic order between the 
two bridges, and none attached to j, or vice versa. The relative portion of the original 
diagram must be a triagulation. 

• No vertices nor edges in _D \ D" exist besides the ones described above. 

Cfr. figure [7]for some examples. Pairs of this form have order 2(n — 3 — \V'\), if 1 < \V'\ < n. 
Remark in particular how this formula is in agreement with the special case V = V{D) and 
D a triangulation, where \E'\ — \E{D)\ ~ 3n — 6. 

A simple proof of Proposition 111.11 is as follows. Of course we have exactly n ~ h > 
vertices in V{D) \ V' . As D^s are connected, each of these vertices have degree at least 1, and 
there exists at least one of these vertices (say, i), adjacent to at least one vertex of V (say, j). 



34 



Figure 7: Two pairs (D, V') giving the leading contribution in the set at n and 
\V'\ fixed. Vertices in V' and not in V' are denoted respectively by black and white 
bullets, and edges m the auxiliary graph D" are in bold. 



In analogy with the drawings in figure [7l we call white and black the vertices respectively in 
V{D) \ V' and in V\ so that there must exist some i white, adjacent to a black vertex j. 

If i has degree larger than 1, then we will see that D cannot have maximal order, because 
we can build a locally modified graph D' which improves the order while remaining in the 
proper (n, h) class. Indeed, there is no loop in D adjacent to j and surrounding i alone. Build 
D' as follows: remove all edges incident on i except for (ij), add the loop on j surrounding i 
alone (this increases the order), then add a subset of edges between j and the other previous 
neighbours of i, up to make D' connected (this can only further increase the order). So we 
get that all white vertices neighbours to at least one black vertex, have overall degree 1. As a 
consequence, all white vertices are leaves, and the second part of the statement is proven. 

For the first part, now consider a face of the diagram D. It visits a number of vertices, 
in cyclic order, possibly with repetitions. Wc have just seen that in this sequence we cannot 
have two consecutive white elements. So, if the face has 4 sides or more, there must be at 
least two non-adjacent terminations corresponding to black vertices. Then, we can add the 
corresponding edge, which is not (planar) multiple (by definition of face), and is thus allowed, 
and this improves the order. From this we deduce that all faces are triangles. □ 

Here we see how the conjecture <"Z„ ^ (t/e)"(A(e" + •••)+ B{e'^ + • • • ) Ine), deduced 
from the analysis of sections [7] and [SJ receives new elements. The case V — % would give a 
contribution of the form A(e"-|-- • • ), while the case \ V'\ = 1, optimal for n > 4 (and marginally 
optimal for n — 4), would give a relative factor g-("-3-|v |) _ g-(«-4)^ compatible with the 
non-analytic contribution i?(e^ -|- • • • ) Ine, and in accord with the fact that this contribution is 
related to the thermodynamic limit, as a non-empty set V' produces an integrand which is not 
blind to the e ^ singularities. At this stage of dimensional analysis, the more subtle possible 
presence of logarithmic factors still does not emerge. 

12 Resummation of Pansy Diagrams 

The defining characteristics of the dominant diagrams get simplified in the case V consists of a 
single "central" vertex. All other n—1 vertices have degree 1, and are connected to the central 
vertex by a single bridge. Then, we have 2n — 5 non-contractible arcs with both terminations 
on the central vertex, the maximal allowed number under the constraint that there are no 
consecutive multiple edges, and producing a triangulation. An example is on the left of figure 
[71 We call these diagrams "pansy" diagrams^ 

Note in particular that, as the adjacency matrix of these graphs is the same for a given 
order n, the integrals T{D, V) are all equal, so we just need to count the diagrams at order n, 

*Both for the clear resemblance with the example in figure, and the fact that a triangulation on the Ricmann 
sphere of our kind is "trilobate" , and some pansy species have petals collected in three main directions. 
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Figure 8: On the left, a typtcal pansy diagram with n = 8. In the middle, the draw- 
ing obtained by inverting the coordinates, and compactifying the radial coordinate. 
On the right, the modification of this drawing which makes clear the bijection with 
triangulations of a polygon with n — 1 sides. 



and evaluate a single integral. 

If one inverts the diagram (in the sense of complex-coordinate inversion z — > for a 
drawing in which the central vertex is at the origin), one understands that pansy diagrams are 
in bijection with the triangulations of regular polygons, an enumeration problem again solved 
by Catalan numbers (cfr. figure [8|). More precisely, there are Cn-s triangulations of a polygon 
with n — 1 vertices (contributing to n-forests), and a cyclic symmetry factor l/(n — 1) should 
be included. So at order n, from the counting of the diagrams we have a factor 



(2n-6)! 



(12.1) 



(n- l)!(n-3)! ' 
and the integral is 

/S jn /■2-K—5 la n—1 

As we know that all contributions for {0} C ]/' C V{D) are subleading, we can freely include 
them, and get 

i{D,m^ E ^(An- r^n7''^ n/^(^^)/(^o'^o)'""'n/(^o,e.). 

{0}<ZV'<ZV(D) •'-^ j=l j j=l 

(12.3) 

We recognize that n—1 integrations are all equal, and factorized, so we can write 

^-^,[e)f{e,ef-^^j^ —^{e')f{e,e')j . (12.4) 

The integral in parenthesis, for /i = 1, is related to the one described in equation (|10.46p . with 

.ie,e)^z-2^l(e + ^^e^y (12.5) 

In particular, for h = 1, at leading order it can be replaced by the numerical constant Ki = 
3 — ''"g^ , corresponding to the limit in which both x and e vanish. More generally, we expect 
that this limit is finite for any value of h, and call it Kh 
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which is the hmit for w — > 2 and g ^ gc oi the quantity (|10.42p in the case I — \. 
The combination x above is useful, as it also appears in f{9, 9) 

f[e,e) = ^{l + 0{e,5)), (12.7) 

so it is convenient to express also ^{6) in these terms 

m = -^x - -^e + Oie^eS, S^)) . (12.8) 



More precisely, it has the structure 

ha 



ith 



m = {l-a)e+-^x + 0{e\ eS, 6')) , (12.9) 



a ^3 + ^, (12.10) 



(12.11) 



and our relevant expression is 

J{D, {0}) = 2--+^A-r^ J[ g (^x + (1 - a) x 
up to subleading terms. So we need to consider integrals of the form 

^,^.^(e,e)-™-^ (12.12) 
for m > —1. Some scale factors overall are easily extracted. Defining 5 ~ \ h(h+i) ^ have 



The expression for Z„, collecting also the combinatorial factor (|12.ip . thus reads 

t-Z„^r (^"7^\ T(A{0}) 
(n — l)!(n — 3)! 

(2n-6)! 1 ^_2n+5r^«-i fha. 



TT (n - 1)! 

For m — —1 and m = (corresponding respectively to = 3 and n — A) the integral Wm is 
not convergent in the limit 5 — > oo, and contains a combination of logarithms that we shall 
discuss: 

5 1 

dT(e + r2)^ =5Ve + 52 + eln((5+ Ve + J2)_ ing. (12.15) 
5 2 

/ dT(e + r2)^^ = 21n(<5+ VTT^) -Ine. (12.16) 

For m > the integral is convergent. The general formula in the limit i5 — > (X) is 

dr(. + r^)— ^ = e- ^"^"^~ \y . (12.17) 
, (2m — 1)!! 
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Our dimensional analysis in powers of e and S is valid for a wide range of values for S (it suffices 
that 5 <C 1, and 6 > Ve), so every term in the result which depends on (5 in a way which is 
not compatible with this arbitrariness must be intended as coming from the analytic part of 
the integral, being this a leading term, as in the cases n = 3 and n = 4, or a subleading one, 
as in the case n > 5 with finite S (as it is legitimate to take 6/y/e ^ 1). Thus, at the aim of 
understanding the non-analytic contribution, these terms can be dropped out, and we have 



( 2m 1)!! ^ m < 

{-my. - 



(m-1)! 
(2m-l)!! 



m > 1 



(12.18) 



Substituting (|12.18[) in (|12.14p . again in agreement with our conjecture, we obtain the expres- 
sions for n — 3 and n — A, that we report here together with the ones (for h odd), already 
derived in equations ()7.17p and (|10.37p 



fZ2{gce-') 



{0{e) analytic, 0{e^)) + 



(0(6^) analytic, 0{e^)) 



2Vh + T 4, 

r-; — e hie 

37r/i4 



Vh+lKh 



' \o{e^) analytic, ©(e^)) + 



167^^/7^+T 



e^lne 
e" Ine 



(0(6^) analytic, ©(e^)) 



e^lne 



(12.19) 
(12.20) 
(12.21) 
(12.22) 



For n > 5, as the overall power of e is negative, the pansy diagrams give a non-analytic 
quantity (in a neighbourhood of e = 0) already without a logarithmic factor, which indeed 
does not occur. 

Substituting (|12.18p in (|12.14p . we get a series for the terms with n > 5 



^rZ„(5,e-^) 



n>5 



n>5 



t_hKh 
e 2 



(n-4)! a(n-5)! 
{n-l)\ ^ 2 (n- 1)! 



(12.23) 



The two series are trivial, and have radius of convergence 1. The behaviour near to this point 
is deducible from the exact expressions 



(n-l) 2 12^ ' 



Ed--: 

ji>5 

^ (1 _ (lizi): = ia;3 ^ _ ^)(2 -7x + Ux') . 



n>5 



(n-l)! 6 
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(12.24) 
(12.25) 



Remark in particular the leading non-analytic behaviour in the full sum (|12.23p . for e \ 



thKh 



n>5 



1 - 



t hK,, 



e 2 



In 1 - 



t hKh 



e 2 



(12.26) 



regardless of the precise expression (|12.10p for a{h). We can interpret the result of this cal- 
culation as describing the curve of critical values gdt] h) in the theory of Spanning Forests, in 
a neighbourhood of i = 0, when the limit (|y(G)|) — > oo has been taken before K{F) oo. 
We already know from elementary means the formula for gc{0;h), equation (|7.5p . The result 
above states that, in the limit described above. 



— In gc{t;h) 



hKh 



(12.27) 
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with Kh a numerical constant, obtained from a single (non-singular) one-dimensional contour 
integral, in equation (jl2.6p . and known for h = 1. 

The whole analysis of the pansy diagrams, besides being a remarkable exact result, shows 
a "topological" fact of this peculiar large-volume limit, namely that, for t, e small and 

e\e,(0:=^, (12.28) 

the partition sum is dominated by forests whose adjacency diagram is compatible with the 
presence of a single gigantic tree, and many small trees. Although this statement is not 
quantitative, it has a clear topological reformulation: at order K{F) — n fixed and in the limit 
(|y(G)|) oo, almost surely there is a single tree T which is neighbour of any other tree T', 
and, for each T', there are edges in G \ F with both endpoints on T, such that the edge, 
together with the unique path connecting the endpoints on T, makes a cycle which encircles 
T'. At the same time, no pairs T' , T" of trees distinct from T are adjacent. 



13 Conclusions and perspectives 

The statistical mechanics of spanning forests on various graphs has two main kinds of criticality, 
in the "probabilistic regime" of t real non-negative, besides the trivial "high temperature" fixed 
point t — > +00. The point t = Q corresponds to a massless theory of a scalar fermion, at the 
light of Kirchhoff Matrix- Tree theorem. Furthermore, at some critical value t* , a percolation 
transition may occur, i.e. for values t < t* there exist trees which occupy a fraction of order 1 of 
the volume {gigantic compents), while for t > t* all trees have a characteristic size, depending 
on t alone and not scaling with the volume. This is the specialization to forests of a feature 
holding more generally for the probabilistic sector of the Random Cluster (Potts) model. It 
is shown numerically in three, four and five dimensions [41], and analytically in the "infinite- 
dimensional" limit of fully connected graphs [42] that t* is finite and non-zero in these cases, 
while it is expected that the arising asymptotic freedom for the model of spanning forests in 
two Euclidean dimensions is due to the fact that the "Kirchhoff" criticality and the percolation 
criticality do coincide exactly for d = 2. 

Indeed, it was also at the aim of understanding rigorously this set of conjectures, that we 
performed the study of the model on Random Planar Graphs, with the aim of combining the 
results with KPZ tools. 

We have however to face a problem in the interpretation of the results of the previous section. 
As we said, we have performed a double limit of (|y(G)|) — > oo and K{F) oo, where the 
first one has been performed before the second one. On the contrary, at least on Euclidean 
lattices, for any finite t, we expect a macroscopic number of components, i.e. \ v(G) \ ) ~ ^i^)^ 
and, in order to have a better understanding on the behaviour of the system, we would like to 
control the double limit (|T^(G)|) , K{F) oo with arbitrary scaling. 

We have seen how, in the Feynman expansion, the accessory parameter V' has a deeper 
meaning than it was legitimate to expect: when, for a given diagram, V' is chosen in order 
to maximize the contribution, we have that vertices respectively in V' and not, correspond 
to gigantic and small trees. So we expect that a control on the microcanonical ensemble 
for n \V'\ would help at the aim of understanding the various limits. This is in a sense 
generalizing the approach of the previous section, where we stated that h = must give an 
analytic contribution, coming from graphs of small size, and we analyzed exactly the leading 
contribution to n = 1 for g ^ gc- 

Proposition 111.11 alreadv identifies the class of leading diagrams at both n and n fixed. We 
plan in the next future to attack the problem of re-summation of diagrams in n at further 
values of n. 
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A Generating functions for /c-trees and Hypergeometric 
functions 

We recall the definitions 

. {{h + l)n)! . / X n A 

n\ihn + l)l 

_ {{h + l)n)\ - Vl^"4' 

^''^"-n!(/in + 2)!' - a; ^,,„ . 

^ ^ n>l 

From the ratio of two consecutive summands 



and 



we have 



Sih 



and 



(A.l) 



(A.2) 



and the definition of generahzed Hypergeometric fmiction [43] 

oo 

pF,(a;b;w) ^a„cc;"; (A.3) 

n=0 

«»+! ^ {n + ai)---{n + ap) ^ 
an {n + bi) ■ ■ ■ {n + bg){n + I) ' 

ao = l; (A.5) 



Ah{uj) = h+iFhiah;hh; Chuj) ; (A. 6a) 



= i {h+iFhiah:h'h;chCj) - 1) ; (A.7a) 
b'. = (f,...,V.^.^)- (A.7b) 

The parameters of functions Ah and have all but one entry in common, the last one differing 
by 1, so they are contiguous as hypergeometric functions (cfr. [43], par. 2.2.1), and this accounts 
for the simple relation p.lSp . 

The hypergeometric function corresponding to Afiito) has been already studied by M.L. 
Glasseil^. In particular, the relation (25) in his communication coincides with our parametric 
solution p.l4p . after that ^/^(cj) has been identified with the hypergeometric function in (|A.6p . 

The hypergeometric functions pFg with p = q + 1 and non-integer parameters do not have 
poles and essential singularities. They have a branch-cut discontinuity between z* — 1 and 
the point at infinity. Accounting for the rescaling constant Ch, we have w* = /{h 1)''+^ — 
2^gc{h) (cfr. definition (|7.5p ). and in particular u* — 1/A ior h — 1, a,s clear from the explicit 
expressions p.Sp and (|3.17p . 

^M.L. Glasser, pers. comm. to the Wolfram Mathworld community, Sept. 26, 2003 
|http : //mathworld. wolfra m . com/HypergeometricFunctioii . html cq. 25. 
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B Comparison with random spanning trees 

The results above for the ensemble of random planar lattice can be compared with small effort 
with the case of random lattices, regardless to the genus. In this case, the only difference is that 
the Catalan numbers, deriving from the combinatorics of planar matchings, must be replaced 
with the number of arbitrary matchings of 2n points, which are (2n— 1)!!. Thus, provided that 
hV is even, we have 

Z,{g)^Y.9''^kvihV+l)ll. (B.l) 

V 

Again we treat separately the two cases of h odd or even. In the case of h odd we have 

M9)-2_,9 ^2V)\{hV + l)l2'-y+^ ' ^ ' 

from which we have the asymptotics 

while in the case of h even we have the formula 

from which we have the asymptotics 

^i(5)-E(i'^^)!(5|^7^)V^ (B.5) 

Remark the expected "entropic catastrophe" , due to the super-exponential number of random 
lattices. As a consequence, the related hypergeometric function is a qFp, with q — p > 1, 
contrarily to what happens in the planar case, in which one deals with p+iFp functions, which 
have a finite radius of convergence. 
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